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Intermediate energy scale physics plays a very important role in non-eqnilibrium dynamics of 
quasi-low dimensional cold atom systems. In this article we obtain the universal scaling relations 
for the generalized reflection coefficient, i.e., the fraction of atoms scattered out of the initial state, at 
intermediate energy scales, scales larger than the trap frequency but much smaller than the effective 
range of the potential, for one and two dimensional harmonically confined geometries. Whenever the 
energy of the cold atoms is commensurate with a transverse energy level, it is shown that the system 
becomes non-interacting. When the difference between the energy of the cold atoms and the given 
transverse energy level, SE, is small compared to the trap frequency, uj, i.e. when S£ = 5E/uj <C 1, 
the reflection coefficient has the universal scaling form R « C^fdE, where C is a constant. The 
power law behaviour and prefactor C appear regardless of the three dimensional scattering length 
and initial conditions. 


I. INTRODUCTION 

The study of cold atoms has flourished due to the vast 
extent of control and tunability available to experimen¬ 
talists. In the laboratory it is now possible to create 
non-equilibrium systems in a controllable manner for any 
dimension, using optical traps, and for any interaction 
strength with Feshbach Resonance mm- Such exper¬ 
iments have led to the study of quench dynamics and 
metastable states EHS], thermalization rates minHii] 
and the dynamics of solitons [H M- These systems 
offer a rich variety of new physics in relation to their 
equilibriated counterparts. One of the major challenges 
of the field of dynamics is the presence of different en¬ 
ergy scales in the problem. Far-from-equilibrium physics 
requires the knowledge of numerous energy scales which 
can differ appreciable, in direct contrast to thermody¬ 
namic physics where only energy scales small compared 
to the temperature dominate. These additional energy 
scales complicate the problems appreciably as the eigen¬ 
states of macroscopic systems at these larger scales are 
extremely complex. 

In the field of dynamics, an interesting issue is whether 
and under what conditions there are universal properties 
at these intermediate energy scales, analogous to the ones 
seen in thermodynamics. Since the number of dynamical 
effects in cold atoms is vast, such a classihcation of dy¬ 
namics into different universality classes would be an ex¬ 
tremely powerful tool. Some efforts in this area have led 
to universal relations for quench dynamics across quan¬ 
tum critical points [SIEHn], turbulence, and soliton dy¬ 
namics in cold atom systems [18) . 

Universal or not, one particularly well studied dynam¬ 
ical effect is the scattering properties of cold atoms in 
confined geometries [iniiniiiMsi- These scattering ex¬ 
periments provide detailed information about the inter¬ 
action between individual atoms, and can be used as a 
starting point for many-body theories. The scattering 
properties of cold atoms are generally examined at low 
energy scales, where the interaction can be parametrized 


by the d-dimensional effective scattering length [21]. In 
the case of harmonically confined geometries, it is pos¬ 
sible to integrate out the transverse degrees of freedom 
and write down an effective low-dimensional model with 
scattering lengths defined in terms of the harmonic trap 
length, a±, and the three dimensional scattering length, 
a uniiinj. Furthermore, this allows one to tune the low¬ 
dimensional interaction strength from zero to infinitely 
positive or negative values by simply varying the three 
dimensional scattering length or the length scale of the 
harmonic confining potential. This effect is known as 
confinement induced resonance m, and is a consequence 
of virtual scattering between different transverse energy 
levels in harmonically confined geometries. This phe¬ 
nomenon has been experimentally verified nmmusi for 
both ID and 2D systems. 

In this article we examine the scattering properties 
of cold atoms and introduce another important conse¬ 
quence of harmonically confined geometries: the univer¬ 
sal properties of scattering dynamics at intermediate en¬ 
ergy scales, or energies larger than the trap frequency uj± 
but much smaller than the effective range of the poten¬ 
tial. In general the scattering dynamics of cold atoms are 
not universal in the sense that they rely on the specihc 
value of the effective low dimensional interaction. How¬ 
ever, whenever the energy of the cold atoms is resonant 
with that of a transverse energy level, the atoms become 
non-interacting. In the proximity of such a resonance, 
physical observables take on universal scaling forms de¬ 
pending only upon the energy difference between the 
atoms and the transverse energy level. The scaling forms 
in this limit are independent of the interaction strength 
and for a wide range of initial conditions. 


II. TWO BODY PROBLEM FOR ARBITRARY 
LOW DIMENSIONS 

Consider two cold atoms in the presence of a confin¬ 
ing potential U{r) = where r_L represents the 


2 


magnitude of the position vector in the transverse di¬ 
rections, interacting via a separable contact interaction 
V^f— r') = gS(f— r'). Such a confinement can be eas¬ 
ily produced in the laboratory mm and acts to reduce 
the original three dimensional system to that of a quasi- 
d dimensional system, where d = 2 — e, and e > 0. In 
the case of harmonic confinement, the center of mass and 
relative motions of these two particles can be separated. 
The interaction will only affect the relative motion of the 
particles, and thus the center of mass motion can be ig¬ 
nored. The Hamiltonian for the relative motion is given 
as: 


H = + U{f) + V{f^, (2.1) 


where m and h have been set to unity. The relative 
Hamiltonian without interactions has a set of eigenstates 
with continuous quantum numbers k for the free d = 2 —e 
dimensions, and a discrete set, n, for the transverse 


energy levels. The energy for a state Iv) 



El, = k'^ + En. The discrete energy levels associated 
with the transverse motion satisfy the harmonic oscilla¬ 
tor Schrodinger equation in the transverse dimensions. 
In this notation, the spectrum of the Hamiltonian is a 
series of bands with quadratic dispersion indexed by n. 

All the scattering dynamics are contained in the T- 
matrix which incorporates the effects of the interaction 
to infinite order. The T-matrix is defined through the 
integral equation: 
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where Gv is the propagator for the non-interacting 
Hamiltonian. Since the interaction potential is separa¬ 
ble, it is possible to solve the scattering problem in its 
entirety. The solution is: 


T{E)-^ = 


r (e /2) 
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(2.3) 

In Eq. (2.3), the free dimensions have been integrated 
out, leaving only the discrete bands a ssoci ated with the 
transverse motion. The sum in Eq. (2.3) is divergent. 
This divergence can be incorporated into a finite renor¬ 
malized coupling constant g. Each band contributes to 
the renormalization process, and it is possible to study 
the flow of the coupling constant as these transverse en¬ 
ergy levels are removed from the system. 

To this end consider a system initially prepared in 
the state |r'o) = no,k\ with energy E. Instead of a 


system with an infinite number of bands, one consid¬ 
ers a system with a finite number of bands with indices 
riQ — s<n<no-l-s and an appropriate renormalized 
coupling constant: 




Figure 1: (Color online) The above figure shows a single 
step in the renormalization procedure. Before the step 
the system is characterized by a finite number of bands 
with indices riQ — s<n<nQ + s and a coupling 
constant g{s), which is shown by the left hand diagram. 
The bands with indices no i s, designated by the thick 
red lines, are then simultaneously integrated out. The 
result is a system with two less bands and a new 
coupling constant g{s — I) as shown on the right. 
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The renormalized coupling constant, g{s), is a function 
of the initial conditions, the energy and, the number of 
bands left in the system, parametrized by s; g(s, E, no) = 
g(s). We require that the T-matrix produces the same 
physics regardless of the number of bands, which implies 
that the T-matrix is invariant for any value of s. This 
invariance condition allows us to calculate the change in 
the coupling constant as s is reduced to zero. The renor¬ 
malization step is defined by integrating out the bands 
uq + s and uq — s (if it exists) by defining a new renor¬ 
malized coupling constant g(s — 1): 


T(E)-^ = 


1 1 

g(s-I) g(s) 
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This procedure is performed until the system only con¬ 
tains one band left, the initial one, and is shown in 
Fig. (0. 

After performing the renormalization step, one finds 
the relation between g{s) and g(s — 1): 
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where I{s) is: 


I{s) = 47ra_L 
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and a_L is the harmonic length of the trap. It is then 
possible to cast this relation into a form that resembles 
the standard renormalization group formulation [13 ng. 
One first defines a dimensionless coupling constant, 


g{s) = y/2uj(s - l)g(s) = -Vs - lg{s), (2.8) 

a_L 

With this definition, in direct analogy to the /3-function, 
we define a{g, s, E), the so called a-function, to describe 
the flow of g: 


^{9,s,E) = 
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In contrast to the /3-function, the a-function depends 
explicitly on the cut-off s. This is due to the fact that the 
harmonic trap explicitly breaks the scale invariance of the 
problem. In spite of this, the a-function still has relevant 
information on the renormalization process in such a sys¬ 
tem. One can show that the only fixed point in Eq. (2.9) 
is the non-interacting fixed point, g{s) = 0, opposed 
to the /3-function which allows a strong coupling fixed 
point in the infra-red limit in the case of one dimension, 
the Tonks-Girardeau gas [29]. This non-interacting fixed 
point can be reached in one of two ways. The first is the 
trivial case where the system is initially non-interacting 
and a = 0 for all steps. However, it is possible to reach 
this fixed point in an atypical manner. 

If the initial dimensionless coupling constant is non¬ 
zero, it is possible to reach the non-interacting case if I{s) 
is divergent for some step s*. At this step the a-function 
has the following form a{g,s*,E) = {s — l)g{s*), which 
implies g{s* — I) = 0, independent of the actual form 
of g{s*) before the step. This universality of the flow is 
shown in Fig.[^(a). The fixed point is thus universal for 
infinite I{s*). If I{s*) is finite, non-universal corrections 


that depend on g{s*) will appear and define the width of 
the universal regime. 

For the specific case of Eq. (2.7), I{s*) diverges 
whenever the energy is commensurate with En» as 
e~^ (dF)”*^^^, where 6E = \E — and n* = uq + s*. 
This creates a sequence of non-interacting universal fixed 
points separated by the band gap in energy space as 
demonstrated Fig. (b) and (c). For large but finite 
/(s*), the dimensionless coupling constant will have the 
universal scaling form: 


~gis* - I) cx 6Eyi^/e. (2.10) 

Corrections to this universal behaviour will appear at 
order SEn* and are depicted in Fig[^(d). In the limit e = 
0, or in two dimensions, /(s*) diverges logarithmically 
and hence the dimensionless coupling constant will also 
scale logarithmically, log“^ |d£|, with corrections of order 
log-2 |(5£-|. 

However, if the kinetic energy of the system is much 
less than the band gap, the initial band and n* will co¬ 
incide. Thus the interaction will renormalize to a finite 
value as the initial band is not included in the renor¬ 
malization procedure. The system then behaves as a 
true low dimensional system which, for one dimensional 
gases, implies the system will form a Tonks-Girardeau 
gas, whereas in two dimensions, the gas will instead 
evolve towards the non-interacting gas. 

This unique structure of the a-funciton has immediate 
implications on the scattering observables of quasi-low 
dimensional gases. This can be deduced by examining 
the fraction of particles scattered out of an initial state, 
i.e. a generalized reflection coefficient. Near these non¬ 
interacting fixed points we can calculate the reflection 
coefficient perturbatively by normalizing the transition 
rate from Fermi’s golden rule [30] with respect to the 
initial particle flux. The matrix element, {f\V\i) can 
be replaced with the renormalized coupling constant as 
it incorporates the effect of all the possible interactions. 
Following the above prescription for the reflection coeffi¬ 
cient one finds: 


R = C (2.11) 

Vc. - no 


where £ is the energy measured in units of the band gap, 
which is 2a; for quasi-one and quasi-two dimensional sys¬ 
tems, and r]{£) is the density of states. The constant C 
is universal and depends specifically on the flow of the 
coupling constant. 

Slightly above a transverse energy level rj{£) diverges 
like with 6£ = £ — n, but is otherwise finite. 

When the energy is commensurate with a harmonic os¬ 
cillator level, Eq. (2.11) will exhibit a universal scaling 
form. For S£ > 0, a partial cancellation occurs between 
the renormalized coupling constant and the density of 
states leading to the universal scaling relation for R: 
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Figure 2: (Color online) The discrete flow of the dimensionless coupling constant, g, is shown above under three 
different conditions: a) fixed initial band no, and total energy £, with varying interaction strength a b) fixed no and 
a with varying £, and c) fixed £ and a with varying uq. In a) the system is prepared with no = 10, and £ « 13. The 
dimensionless coupling constant will universally flow into the non-interacting hxed point when s = 3 regardless of a. 
A similar situation is shown in b). For the same no and a, the initial energy is varied in units of the band gap. g 
once again universally renormalizes to zero at the step s with £ = uq + s. Finally, in plot c) a and the energy £ are 
fixed with £ = 12. The flow is shown for three different choices of the initial band index uq. g then renormalizes to 
zero for all n* ^ tiq. Only in the case n* = no will g renormalize to some non-zero value dictated by the true low 
dimensional system. In part d) the value of the coupling constant after renormalization | 5 (s* — 1)| is shown for 
different values of d£ with no = 10, s = 2, and £ = 12 + 5£ for a quasi-one dimensional system. For finite 5£, there 
are corrections to the universal form of |g(s* — 1)| oc of order 5£. 


R = CSr/'^le. (2.12) 

When 5£ < 0, the density of states is finite and as a 
result there is a different scaling relation: 

R = CSr/e. (2.13) 

Finally, when e = 0. The case of a quasi-two dimensional 
systems, the density of states is a constant, and there is 
only one scaling relation for R\ 

R = C\og-‘^\5£\, (2.14) 



Scaling Exponent 


< 0 

SS >0 

ID 

0.5 

1 

2D 

2 

2 


Table I: Universal scaling exponents for quasi-one and 
quasi-two dimensional systems when the energy, in 
units of the transverse energy scale, £, is commensurate 
with a transverse energy level. In the case of ID 
systems, the scattering observables scale as \5£\^, for 
some exponent x. Similarly for 2D, the observables 
scale as log~“ |i5£|. The scaling coefficients are shown 
for the two possibilities; when the energy approaches 
some excited band from below or above. 


valid whenever \5£\ 1. The scaling exponents for these 

relations are found in Tablej^for quasi-one and quasi-two 
dimensional systems. 

We now compare Eq. ( |2.11[ ) with the exact results ob¬ 
tained for quasi-one and quasi-two dimensional systems. 
For quasi-one dimensional gases, a calculation of the re¬ 
flection coefhcient for arbitrary energies is given by m- 


(2.15) 

where C(l/2,a:), is the Hurwitz zeta function |31j . 

For positive £, the Hurwitz zeta function diverges peri¬ 
odically whenever £ is equal to an integer, or equivalently, 


R = il{£) ^+ai/2,-£) 
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when the energy is near to a harmonic oscillator level, as 
Upon expanding this exact solution around the 
non-interacting fixed points one finds the same universal 
behaviour and scaling exponents as the renormalization 
analysis. 

Our perturbative analysis of the reflection coefficient 
holds for every non-interacting fixed point. When n* = 
riQ the system will renormalize into a strong coupling 
regime. In this case the exact solution gives: 


1-R'. 




C(i/2)), 


(2.16) 


fully consistent with the predictions from the a-function. 

For a quasi-two-dimensional system, it is the cross 
section cr that exhibits periodic divergences whenever 
£ is an integer. The exact solution for a is found in 
Ref. [TU]. This solution has the same scaling behaviour 
as found by Eq. (2.111. The flow for quasi-two and two di¬ 


mensional systems are both described by non-interacting 
fixed points. 


III. CONCLUSION 

In this analysis we considered a discrete renormaliza¬ 
tion procedure to study the contact interaction between 
two cold atoms in a confined low dimensional geome¬ 
try. The renormalization procedure was performed by 
integrating out each band describing a transverse energy 
level. If the energy of the system was commensurate 


with the bottom of a new band {n* ^ ng), the interac¬ 
tion then renormalized to zero, resulting in total trans¬ 
mission and a non-interacting theory. Near these critical 
points, observables take on universal scaling forms and 
their scaling exponents were determined. This approach 
is valid for a wide range of initial conditions, and as a 
result is quite ubiquitous. The only caveat is when the 
kinetic energy of the initial system is much less than the 
band gap 2oj, i.e n* = uq. In this case the interaction 
renormalizes to a constant, and the physics is dominated 
by the true low dimensional system. For quasi-one di¬ 
mensional systems, this observation leads to the strong 
coupling Tonks-Girardeau gas [2^ . 

The question of universal dynamics is still in its in¬ 
fancy, but has resulted in several fascinating results in 
quench dynamics HUMI] and cold atom systems |I8) . 
It remains to be seen whether other dynamical systems 
can fall into universality classes. The number of different 
dynamical effects is quite large, and the conditions un¬ 
der which universal dynamics can be obtained is a per¬ 
tinent question. If an understanding of when universal 
dynamical effects are present, the categorization of dif¬ 
ferent universality classes can be conducted thoroughly 
and explicitly tested in the laboratory. 

The authors would like to thank Shao-Jian Jiang for 
useful discussions. This work was funded by the Cana¬ 
dian Institute for Advanced Research (CIFAR) and the 
Natural Sciences and Engineering Research Council of 
Canada (NSERC). 

In preparing this manuscript, a similar work was done 
in Ref. (35] using numerical solutions of an exact K- 
matrix formalism. 
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